Simple model for decay of superdeformed nuclei 
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Recent theoretical investigations of the decay mechanism out of a superdeformed nuclear band have 
yielded qualitatively different results, depending on the relative values of the relevant decay widths. 
We present a simple two-level model for the dynamics of the tunneling between the superdeformed 
and normal-deformed bands, which treats decay and tunneling processes on an equal footing. The 
previous theoretical results are shown to correspond to coherent and incoherent limits of the full 
tunneling dynamics. Our model accounts for experimental data in both the A ~ 150 mass region, 
where the tunneling dynamics is coherent, and in the A ~ 190 mass region, where the tunneling 
dynamics is incoherent. 
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One of the most interesting recent discoveries in 
nuclear-structure physics is the existence of superdefor- 
mation for nuclei in the mass A ~ 150 and A ~ 190 re- 
gions. So far, a consistent theory regarding the decay out 
of a superdeformed (SD) rotational band into a normal- 
deformed (ND) band has not been achieved. Most of the 
early work |l|-|6|] on this problem attributed the decay-out 
process to a mixing of the SD states with ND states of 
equal spin. Decay out of the SD band sets in at a spin 
Io for which penetration through the barrier between the 
SD minimum and the ND minimum is competitive with 
the E2 decay within the SD band. A statistical model 
was used [||!| to describe the ND states, and the de- 
cay out of the SD band was determined as a function 
of the decay widths Ys and IV in the SD and ND po- 
tential wells, respectively; the spreading (or tunneling) 
width r^ through the barrier; and the average spacing 
Dm of the ND states. Under the assumption that the 
ND states form a continuum on the scale of the other 
energies in the problem, the spreading width was found 
using Fermi's golden rule to be pj 



Y l = 2?r(V 2 )/D 



n- 



(1) 



where (V 2 ) is the mean square of the coupling matrix 
elements V a p connecting the SD and ND states. T^ mea- 
sures the strength of the coupling between the SD and 
ND states. In Refs. J2|||, it was assumed that IV < Is 
and that Y^- /Dm ~ 1, i.e., that the coupling between the 
SD and ND states is relatively strong. 

Quite recently, a different approach to this problem 
has been reported 0. In this approach, the reduction 
factor Fs of the intraband transition intensity is calcu- 
lated directly as a function of the spreading width T^ 
and of the intraband E2 width Yg. Their final result for 
Fs is shown to be independent of the statistical El decay 
widths IV of the ND states, provided that IV 3> T*, Yg. 
Since the publication of the later result, it has been diffi- 
cult to reconcile the predictions of these two calculations, 
because their final results do not depend upon the same 
parameters. 



The purpose of the present manuscript is to formu- 
late a simple two-level model for this problem, so as to 
study in detail the dependence of the decay-out process 
on Ys, IV, and I>. It will be shown that the results of 
both previous investigations of the decay-out process can 
be obtained in certain limits, depending on the relative 
sizes of these widths. It will also be shown that the ap- 
propriate expression for the spreading width is not Eq. 
(0), but 



T i = 



2YV 1 



S2 ' 



(2) 



where f = (Yg + IV)/2, V is the matrix element con- 
necting the SD state of interest with the single ND state 
with which it mixes most strongly, and A is the energy 
difference of these two states. It will be shown that Eq. 
(jy) is indeed the correct mean value of Y- 1 over a statisti- 
cal ensemble of nuclei in the limit V <C Y, in agreement 
with Ref. M. However, the fluctuations in Y^ are typi- 
cally much larger than its mean, indicating that Eq. (0) 
must be used to describe the decay out of a particular 
superdeformed state. 

Motivated by the experimental fact M that IV, Ys <C 
Dm in the A ~ 190 region, we consider an effective two- 
level system consisting of the superdeformed state S and 
the normal-deformed state N to which it couples most 
strongly. The Hamiltonian of the system is 



H = H + H L 



(3) 



where Ho describes the electromagnetic decay processes 
IV and Ys within the ND and SD bands, and 



#o = £sclc s + e N c N c N + V (c s c N 



c N c s 



(4) 



describes the effective two-level system, including tunnel- 
ing through the barrier separating the SD and ND states. 
Here c s and c N are creation operators for the superde- 
formed state S, of energy £5, and the normal-deformed 



state N, of energy en, respectively. Without loss of gen- 
erality, the tunneling matrix element V is chosen posi- 
tive via an appropriate choice of the relative phase of the 
states S and N. 

In order to include both the coherent "Rabi oscilla- 
tions" due to V and the irreversible decays T$ and IV, 
it is useful to consider the retarded Green's function 



G ij (t) = -ie(t){{c i (t),c](0)}) 
and its Fourier transform 



/oo 
dtGij(ty 
-oo 



iEt 



(5) 



(6) 



The Green's function of the tunneling Hamiltonian Hq 
satisfies 

G Q 1 (E) = l(E + iQ + )-H Q , (7) 

where 1 is the unit matrix. In the \S), \N) basis, one has 



Go 1 ^) 



E-e s + iO + -V 

-V E - e N + i0+ 



(8) 



The full Green's function, including decay processes, may 
be calculated from Dyson's equation 



G- 1 = Go 1 - S, 



(9) 



where £ is the self-energy matrix describing the decay 
processes T$ and IV induced by Hjj. The simplest 
ansatz for £ is M 



£jVS ^NN 



-iTs/2 

-iT N /2 



(10) 



Using Eqs. (ph and (JIG), one can solve Dyson's equation 
to obtain the full retarded Green's function of the two- 
level system, 

q = ( Gss Gsn 
\ Gns Gnn 

= [(E -e s + iTs/2) {E-e N + iT N /2) - V 2 ' 

E - e N + iT N /2 -V 



-V 



E-E S + iTs/2 



(11) 



From G, all information about the dynamics of the sys- 
tem and the branching ratios of the decay processes can 
be obtained. 

Let us first study the dynamics of the coupled SD-ND 
system. Assuming the nucleus starts out at time zero in 
the superdeformed state \S), the probability that the nu- 
cleus is in state \S) at a later time t is Ps(t) = \Gss(t)\ 2 - 
The probability that the nucleus is in the normal state 
\N) at time t is P N (t) = \G NS (t)\ 2 . P s (t) and P N {t) may 
be calculated straightforwardly from the Fourier trans- 
form of Eq. (pd|). The general result for P/v(t) is 



P N (t) = 



2V 2 



-Ti 



(coshojit — cosoj r t) 



(12) 



where 



= v/4U 2 + (A-zr') 2 , 



(13) 



with A = e N - e s and V = (rjv - r s )/2. The general 
expression for Ps(t) is rather lengthy. 

The tunneling dynamics is particularly interesting 
when the energy difference A = 0. There are then two 
qualitatively different dynamical regimes, depending on 
the relative size of the tunneling matrix element V and 
the difference in decay rates T'. For 2V > \T'\, the tun- 
neling dynamics is coherent, and one finds 



Ps(t) 

P N (t) 



er Tt 
W 2 



cos 



LV Q t V 



^12 



H sinwoi + 



U) t 



—Vt 

e sin 



2 ujq 
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uJ. 
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where the Rabi frequency loq is 



uj q = Uv 2 -T 



/2|l/2 



(14) 



(15) 



For 2V < |r'|, on the other hand, the tunneling dynamics 
is incoherent, and 



Ps(t) - 



-n 



Wo 



2 uj t , r' r'^ . 2 w oi 

cosh 1 sinhwot H r smh 



u,', 



II 



4^2 



P N (t) = ^ e -"sinh 2 ^. 



w, 



(16) 
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For r" = and A^O, the tunneling dynamics is coher- 
ent, given by Eq. @ with uj -> (AV 2 + A 2 ) 1 / 2 . For 
r' / and A / 0, the tunneling dynamics has both 
coherent and incoherent components [c.f. Eq. (fL2|)], the 
coherent component being suppressed for large V and/or 
large A. 

The dynamics of the system is similar to that of 
the two-level system with dissipation, investigated by 
Leggett and collaborators M; the principal difference is 
that we consider a two-level system in which the total 
number of particles is itself time-dependent. The physical 
origin of the imaginary self energy £ is virtual transitions 
of the nucleus to lower-lying states and back again, which 
alter the state of the electromagnetic environment. This 
is analogous to the coupling of the two-level system to a 
bath of bosonic excitations considered in Ref . |1 . If the 
environment couples with equal strength to the states 
S and N, i.e., if T$ = r^v = f , the Green's function 
factorizes quite generally |Hj: G(t) — e~ rt / 2 Go{t), and 
the nucleus undergoes Rabi oscillations with frequency 
lu = (AV 2 + A 2 ) 1 / 2 between the states S and TV. The 
nucleus is in a coherent superposition of states, which 
decays at an overall rate f to lower-lying states. How- 
ever, if the environment couples with different strengths 
to the states S and N, i.e., if f ^ 0, coherent tunneling 
between S and N is suppressed since the environment 
"measures" which state the system is in. For A = and 
< |r'| < 2V, the dynamics described by Eq. (|3) is 



qualitatively similar to the case T' = 0, but the Rabi fre- 
quency is reduced to the value given in Eq. (|l5|). If the 
difference in coupling exceeds the critical value |r'| > 2V 
for A = 0, the coherent superposition of S and N is de- 
stroyed altogether, and the dynamics is overdamped. As 
in the model of Ref. || , there are both coherent and in- 
coherent components of the dynamics when both V =/= 
and A^O. 

Let us next turn our attention to the decay branching 
ratio, which is the experimentally measurable quantity. 
When the nucleus is in the state S, it decays at a rate Tg 
to a lower superdeformed state, and when it is in state 
N it decays at a rate Tjv to a lower-energy state in the 
normal-deformed band. Thus, the time-dependent rates 
to decay in the S and N channels are 

f s (t) = r s p s (t), f N (t) = r N p N (t). (17) 

The fraction Fn of nuclei that decay via El processes in 
the normal-deformed band is just pT| 



that Ts is an E2 decay and T^ is an El decay. In the 
limit T N > Ts, Eq. fll9| ) simplifies to 



lim 



Fff — 



r' 



r s + ri' 



(23) 



where T^ is given by Eq. (g). Eq. ( j23] ) is identical to 
the principal result of Weidenmullcr, von Brcntano, and 
Barrett M, although our expression for T^ differs from 
that of Ref. [[7j. Note that, in contrast to the argument of 
Ref. M , no assumption has been made about the relative 
size of V and f in deriving Eq. ( |23|) from Eq. (|lj) . How- 
ever, the interpretation of I^ as a tunneling rate is only 
justified when r^/f < 1; for larger values of V, the tun- 
neling dynamics described by Eq. (n3) is more complex, 
though the integrated rate still obeys Eq. (23), provided 
Tat » T s . 

Eq. (g) for r^ is also the expression one would obtain 
from a correct application of Fermi's golden rule in the 
limit V < f < D N : 



Fn = 



f °°dtf N (t) 



J oc dt[T N (t) + r s (t)} 



r°° re 

= T N / dtP N {t). (18) r i =2nV 2 / 

Jo J- t 



dEp s (E) PN {E), 



(24) 



This integral may be evaluated to obtain the central re- 
sult of this paper, 



Fn — 



(1 + T N /T S )V 2 



A 2 +T 2 (1 + W 2 /T N T S )' 



(19) 



The fraction of nuclei decaying via E2 processes within 
the superdeformed band is Fs = 1 — Fn ■ In Ref. JjJ , Fs 
was denoted by F. 

Let us next consider some limits of Eq. (|T^) . In the 
limit of very strong coupling of the states S and N, V 3> 
r, one finds 



,. p, Ljy(rs +Tjv) 

v/ ™oo N ~ (T S + r N ) 2 + T S T N (A/V)2 ' 



(20) 



This is equivalent to the result of Vigezzi, Broglia, and 
Dossing B,p| for the case where only a single SD and ND 
state mix: 



T^vbd 
*N 



2-^, (i _ |„ I2^r., j_ u I2i 



(i-|c CT | 2 )rv + M 2 iy 



(21; 



where c± — {±\S) are the mixing amplitudes of the eigen- 
states |±) of Ho with \S). From Eq. ([|), we have 



|c±| 2 = [l + (x± Vx 2 + 1) 



(22) 



with x = A/2F. Inserting Eq. (§|) into Eq. @, one 
indeed recovers Eq. (p0|). Thus the result ( pi| ) of Rcfs. 
[g,g| is seen to be a limiting case for V/T —* oo (i.e., for 
fully coherent tunneling) of our general result, Eq. (|l9|). 
Another limit was considered in Ref. , namely Tn ^ 
Ts, TK In this limit, the tunneling dynamics is incoher- 
ent. The assumption Ts <S Tn is motivated by the fact 



where the lifetime-broadened densities-of-states of the SD 
and ND levels are 



Ps(E) = 
Pn(E) = 



T s /2n 



{E-es) 2 + T 2 s /^ 

Tn/2h 
{E-e N y+T 2 N /t 



Evaluating the integral in Eq. (J2J), the expression (g) 
is obtained. The level-spacing Dn in the ND band is 
irrelevant if V <€. Dn, since V only mixes the state S 
and the single state A^ which is closest to it in energy in 
that case, as we have assumed. 

The branching ratio dl^) depends strongly on the en- 
ergy difference A = sn — £s, which in turn depends sen- 
sitively on the microscopic Hamiltonian of the particular 
nucleus under investigation. In order to eliminate this 
parameter dependence, one practice which is employed 
is to calculate the average of Fn over a statistical (GOE) 
ensemble of similar nuclei. In the limit of incoherent tun- 
neling V <§; r,5, Tn, Eq. (|9|) may be integrated over A 
to obtain 



(Fn) 



^f N (a) = n 

D N K ' T S 



(25) 



where (r^) is given by the right-hand-side of Eq. (Q), 
in a gree ment with Ref. H. However, it is clear from 
Eq. ( |19|) that Fn typically deviates significantly from its 
mean value. For instance, the mean square of Fn is much 
larger than the square of (Fn) when Dn 3> f : 



(F 2 N ) D N (V<) 



(Fn)* 2irT(V 2 ) 



2\2 ' 



(26) 



Thus, it would be preferable to compare experimental re- 
sults directly with Eq. ([19|), rather than with its ensemble 
average. 

In Table I, we show some experimental data for nu- 
clei in the A ~ 150 and A ~ 190 mass regions. So far, 
little data is available in the A ~ 150 region, with only 
estimates pjl^l f° r the widths for 152 Dy. From the num- 
bers given in Table I, we observe that the nuclei in the 
A ~ 190 mass region have T+, Ts <C Tjv. The dynamics 
of S — > N tunneling in these nuclei is thus incoherent, 
and the appropriate branching ratio is given by Eq. (23), 
in agreement with Ref. |j. On the other hand, for 15 ^Ty, 
r^ 3> Ts, IV, indicating coherent S ^ N tunneling. The 
measured value 11211 of F$ = 0.51 for 152 Dy(26) is con- 
sistent with Eq. (|2C|), using the values of Ts and Tjv in 
Table I and assuming V/A ~ 1, in accord with the the- 
ory of Refs. |l|-|(|. Our general result (|19|) is consistent 
with the data in both the A ~ 150 and A ~ 190 mass 
regions and unifies these two complementary theoretical 
approaches. 

A recent paper |13) shows that the S — > N tunneling 
rate may be enhanced by several orders of magnitude 
if the ND states are chaotic at the moment of the decay 
out. These results are not inconsistent with our two-level 
model, but would simply imply an enhancement of the 
tunneling matrix element V. 

In conclusion, we have shown that a simple two-level 
model, which treats decay and tunneling processes on an 
equal footing, can explain the apparently disparate previ- 
ous theoretical results, i.e., Refs. |l|-g] versus Ref. [0, for 
the decay out of a superdeformed band. These previous 
results are shown to correspond to the coherent and in- 
coherent limits, respectively, of the tunneling dynamics, 
and are special cases of our general result, Eq. (|l9|). We 
remark that it is straightforward to extend our method 
to treat an SD state coupled to an arbitrary number of 
ND states. 
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TABLE I. Widths and level spacings for a number of nu- 
clei, deduced from the data of Refs. [2,8], following the analy- 
sis of Refs. [2,12] for A ~ 150 and of Ref. [7] for A ~ 190. The 
spin values of the decaying states are given in parentheses. 



Nucleus 


F N (meV) 


D N (eV) 


T s (meV) 


T x (meV) 


iM Dy(26) 


10-20 


3-10 


~ 11 


900-3000 


152 Dy(24) 


10-20 


3-10 


~7.6 


900-3000 


192 Hg(12) 


10.3 


34 


0.116 


0.018 


192 Hg(10) 


10.3 


30 


0.054 


0.544 


194 Hg(12) 


18.1 


92 


0.144 


0.097 


194 Hg(10) 


18.4 


79 


> 0.047 


>0.89 


194 Pb(10) 


1.6 


1699 


0.066 


0.011 


194 Pb(8) 


1.7 


1549 


0.028 


0.009 
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